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j—{ Abstract 
> 

— In this article we explore the ability of dynamical systems tools 

to describe transport in oceanic flows characterized by data sets mea- 
sured from satellite. In particular we have studied the geometrical 
skeleton describing transport in the Kuroshio region. For this pur- 
pose we have computed special hyperbolic trajectories, recognized as 
distinguished hyperbolic trajectories, that act as organizing centres 
of the flow. We have computed their stable and unstable manifolds, 
^ ^ and they reveal that the turnstile mechanism is at work during several 

^ spring months in the year 2003 across the Kuroshio current. We have 



^ found that near the hyperbolic trajectories takes place a filamentous 

transport front-cross the current that mixes waters at both sides. 



1 Introduction 



Transport processes across the ocean have an important impact in global 
climate [Bower et al.(2009)|[Ka"tija fc Dabiri(2009)| , as oceans are important 
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reservoirs of heat and carbon dioxide which are the most important driving 
forces of chmate change. Understanding and describing transport processes 
on the ocean have been a challenging task over the past years, since great 
currents over the ocean such as Gulf stream or Kuroshio that seem smooth, 
river-like streams in their Eulerian description, present a messy pattern on 
the float tracking Bower et al. (2009) , [Lavender et al.(2000)| . Understanding 
such patterns is an essential ingredient for predicting the effects of global 
warming in Earth's future climate. 

Dynamical systems theory has provided in recent years a framework for 
describing transport in time dependent flows. It has been useful for find- 
ing transport routes in complex flows, which means it has successfully found 
order structures from erratic patterns. Typical tools in dynamical systems 
theory are Lyapunov exponents and invariant manifolds. These concepts are 
defined strictly on infinite time systems, however realistic flows, like those 
arising in geophysics or oceanography, are not infinite time systems and for 
their description, finite time versions of the definition of Lyapunov expo- 
nents, such as, finite size Lyapunov expone nts FSLE [Aurell et al.(199"7)] 
and finite time Lyapunov exponents FTLE |Nese(1989) Haller(2001)| are 
used. Invariant manifolds are objects that in the classical theory of dy- 
namical systems are asymptotically defined for special trajectories such as 
fixed points or periodic orbits with hyperbolic stability. Recent articles 
by |Ide et al.(2002)j|Ju et al.(2003)|[M^rid fc Mancho(2009)| have provided 
definitions that generalise these reference orbits for time dependent aperiodic 
flows. These special trajectories are called distinguished trajectories. The 
definition given in [Madrid fc Mancho(2009]] is proven to encompass fixed 
points and periodic orbits as particular cases. Separation points (SP) on 
boundaries such as, detachment and reattachment points are also particular 
examples of distinguished hyperbolic trajectories (DHT). The computation of 
SP has been also addressed in [Haller(2004)[ however the definition provided 
Madrid fc Mancho(2009)] provides a systematic way for finding these spe- 



m 



cial trajectories not only on boundaries but also in the interior of the flow. 
DHT and their stable and unstable manifolds are the essential dynamical 
tools we will apply in this article for describing transport in geophysical 
flows. 

Lagrangian and Eulerian observational datasets have held great promise 
of providing deep insight on oceanic transport. Drifters observations have re- 
vealed fluid exchange Bower &; Rossby(1989) Lozier et al.(2003j] in the Gulf 
Stream. An important step forwards in the description of exchange of trajec- 
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tories across the Gulf Stream was the meandering kinematic model proposed 
by |Bower(1991)| . Studies such as Lozier et al.(2003)] have succeeded in de- 



scribing transport in simplified models of the Gulf Stream by means of lobe 
dynamics |Bower(1991)| . Recently, ocean-observing satellites are entering a 
new operational era and naturally occupy an important place in programs de- 
signed to manage and predict ocean and climate change. These instruments 
may provide a near real time processing system of satellite altimeter data. In 
this article we describe cross-frontal transport along a major ocean current, 
the Kuroshio current, based on a Lagrangian description of observed altime- 
ter velocities, over a period of the year 2003. Our description uses dynamical 
objects such as the stable and unstable manifolds of distinguished trajecto- 
ries. We detect the turnstile mechanism across a part of the Kuroshio current. 
The turnstile survives persistently from April 3 to May 26, 2003. This mech- 
anism had already been found to be present in realistic numerically produced 
data sets modeling the North Balearic front Mancho et al.(2008)| . In this ar- 



ticle we discuss similar ideas with the challenge of dealing with altimeter data. 
Lagrangian Coherent Structures (LCS), as the use of finite time Lyapunov 
exponents to compute transport issues in geophysical fiows is referred to, or 
FSLE have also recently applied to describe dynamics in altimeter datasets 
Beron-Vera et al.(2008)[ [Rossi et al.(2008)| [Rypina et al.(2009)| , although 



these studies are focused in areas of the ocean surface different from the 
Kuroshio current. 

The structure of the article is as follows. In Section 2 we describe the 
analysed data and to what extent the velocity field satisfies the divergence 
free approach. Section 3 describes the details of the interpolation used in the 
data and provides the equations of the continuous dynamical system used on 
our analysis. Section 4 describes the computation of the distinguished hyper- 
bolic trajectories and of their stable and unstable manifolds for the system 
built up in Section 3. Section 5 describes the turnstile mechanism governing 
the dynamics across a current that persist on the data amid latitudes 32°N 
and 40°N and longitudes 155°E and 165°E during the months of April to 
May 2003. Finally Section 6 details the conclusions. 

2 The altimeter dataset 

In this work we analyze transport across the Kuroshio current from April 3 
to May 26, 2003, spanning a time interval of 54 days. 
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We use daily maps of surface currents computed at CLS in the framework 
of the SURCOUF project Larnicol et aL(2006)| from a combination of ahi- 
metric sea surface heights and windstress datain a two-steps procedure: On 
one hand, multimission (ERS-ENVISAT, TOPEX- JASON) altimetric maps 
of sea level anomaly (SLA) are added to the RIO05 global Mean Dynamic To- 
pography Rio & Hernandez (2004) Rio et al.(2005)| to obtain global maps 
of sea surface heights from which surface geostrophic velocities {ug, Vg) are 
obtained by simple derivation. 

gdh 

9 dh 

where / is the coriolis parameter and g the gravitational constant. 

On the other hand, the Ekman component of the ocean surface current 
[uek, Vek) is estimated using a 2-parameter model: = he^^r where h and 9 
are estimated by latitudinal bands from a least square fit between ECMWF 
6-hourly windstress analysis r and an estimate of the Ekman current obtained 
removing the altimetric geostrophic current from the total current measured 
by drifting buoy velocities available from 1993 to 2005 (the method is de- 
scribed in further details in Rio fc Hernandez(2003)| ). 

Both the geostrophic and the Ekman component of the ocean surface 
current are added to obtain estimates of the ocean total surface current that 
are used in the Kuroshio area for the present study. 

Fig. [T] shows a sample of the velocity field on April 4, 2003. We have 
marked with a line the current across which we study transport. The presence 
in the nearby of two counterclockwise eddies suggest the presence of two 
hyperbolic points at both sides of the line that could maintain the turnstile 
mechanism. We will show in Section 5 that this is truly case. The interaction 
of these type of eddies with the main current has been pointed out to trigger 
the Kuroshio meander formation Waseda et al.(2003)| . 

Our work considers that in the area under study the particle motion is 
restricted to the ocean surface, i.e, there is no significative vertical motions. 
The altimeter-derived velocity field is two-dimensional, however since the 
Ekman correction has been introduced, it is not guaranteed that the data is 
divergence free. We have computed the divergence of the velocity field in the 
whole period and area under study and we have verified that despite these 
corrections the divergence free approach is still valid. 
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Figure 1: Velocity field of the Kuroshio current on April 4, 2003. A line 
marks a piece of current across which transport is studied. Maximum values 
of the velocity field are about 3.65 m/s. 

3 The equations of motion 

The equations of motion that describe the horizontal evolution of particle 
trajectories on a sphere is 

# _ u{(j),X,t) 



dt i?cos(A) 
d\ v{(f>,\,t) 
'dt ~ R 



(4) 



where u and v represent the eastward and northward components of the 
altimetry surface velocity field respectively described in the previous sec- 
tion. Particle trajectories must be integrated in equations ([3])-(|4]) and since 
information is provided just in a discrete space-time grid, a first issue to 
deal with is that of interpolation of discrete data sets. A recent paper by 
Mancho et al.(2006)| compares different interpolation techniques in tracking 



particle trajectories. Bicubic spatial interpolation in space [Press et al.(1992)| 



and third order Lagrange polynomials in time are shown to provide a com- 
putationally efficient and accurate method. We use this technique in our 
computations. However we notice that bicubic spatial interpolation in space 
as discussed in Press et al.(1992)] requires an equally spaced grid. Our data 
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input is expressed in spherical coordinates, and the grid is not uniformly 
spaced in the latitude coordinate. In order to interpolate in an uniformly 
spaced grid, we transform our coordinate system (A, 0) to a new one (/x, 0). 
The latitude A is related to the new coordinate fi by 

/i = ln|secA + tanA| (5) 

Our velocity field is now on a uniform grid in the (/i, 0) coordinates. The 
equations of motion in the old variables are transformed to a new expression 
in the new variables, 

d(f) _ u{(l),fj.,t) 



dt -Rcos(A(/i)) 

dt -Rcos(A(/i)) 
where A(/x) is obtained by inverting Eq. (|5]), i.e. 



(7) 



A = ^ - 2atan(e"^) . (8) 

Once trajectories are integrated from these equations, for presentation pur- 
poses one can convert /i values back to latitudes A just by using (Isl). 



4 Distinguished Hyperbolic Trajectories and 
Manifolds 

Distinguished hyperbolic trajectories and their unstable and stable manifolds 
are the dynamical systems objects used to describe and quantify transport. 

Fig. [T] marks with a line the eulerian feature across which we study 
transport. The current, similarly to that studied in Mancho et al.(2008)| 



CouUiette and Wiggins(2001)| , flows eastward. First we identify a DHT in 



the western part of the flow, for which we compute the unstable manifold 
and a DHT in the eastern part, for which we compute the stable manifold. 

Computation of distinguished hyperbolic trajectories for aperiodic flows 
has been discussed in |Ide et al.(2002)| |Ju et al.(2003)| [Mancho et al.(2004)j 



Madrid fc Mancho (2009)]^ The approach taken in this article is that of 



Madrid fc Mancho(2009)] , which is based on the function M defined as fol- 



lows. Let (0(t),/i(t)) denote a trajectory of the system (|6])-([7]). For all initial 
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Figure 2: Structure of the function M for r = 15, a) in the neighbourhood 
of the western DHT, DHTvi/, on March 13, 2003; b) in the neighbourhood of 
the eastern DHT, DHTg, on June 19, 2003. 

conditions (0*, /i*) in an open domain B of the ocean surface, at a given time 
t*, consider the function M {(/)*, fi* ,t*)r : i3 x t — > M 

For an initial condition at t*, the function M measures the length of 

the curve outlined by a trajectory from t* — r to t* + r on the plane (0, /i). 

We discuss in more detail the numerical evaluation of M as defined in Eq. 
Q. Trajectories {(j){t) ^ iJL{t)) of the system (|6])-([7]) are obtained numerically, 
and thus represented by a finite number of points, L. A discrete version of 
Eq. g is: 




where the functions and represent a curve interpolation parametrized 
by p, and the integral 
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is computed numerically. Following the methodology described in [Madrid fc Mancho(2009)] 
we have used the interpolation method used by |Dritschel(1989)| in the con- 
text of contour dynamics. To compute the integral (11) we have used the 



Romberges method (see [Press et al.(1992)| ) of order 2K with K = b. 



The function M is defined over an open set, so it does not necessarily 
attain a minimum, but if it does, the minimum is denoted by min(Mr). In 
Madrid fc Mancho(2009)] it is found that the position of the minimum of M 



at each t* is a function of r. For r >> 0, the minimum converges to a fixed 
value called limit coordinate. Distinguished trajectories are those that for 
a time interval, pass close enough (at a distance e, typically within numer- 
ical accuracy) to a path of limit coordinates. Computing limit coordinates 
becomes a method for finding distinguished trajectories because it is quite 
usual that there exist trajectories passing near these paths. 

We have computed limit coordinates at opposite sides of the line depicted 
in Fig. [TJ A typical structure of the function M for r = 15 appears in 
Fig. [2] The color code reads as follows: maximum M values are dark red, 
while minima are dark blue. As discussed in Mendoza fc Mancho(2009)] 
the important information provided by M are its minima and its singular 
features, so no specific colorbar is required. Fig. |2^) shows the contour plot 
near the western DHT, DHTj^, in March 13, 2003; Fig. |2}d) depicts the 
contour plot of M for the eastern DHT, DHT^;, in June 19, 2003. The 
time evolution of two paths of limit coordinates, DHTyi/ and DHT^, at the 
West are displayed in Fig. |3| Initial conditions on each limit path evolve 
staying near it for a certain time interval. However to prove that a trajectory 
stays close to the path in long time intervals it is a difficult task because 
it is not possible to obtain the whole DHT by direct integration methods. 
As discussed in [Madrid fc Mancho(2009")] , DHTs are elusive and any initial 
condition in their neighborhood, eventually leaves it through the unstable 
manifold. If paths in Fig. [3] are trajectories they should satisfy Eqs. (l6])-([7]). 
However verifying this requires the computation of the time derivative of the 
paths of limit coordinates and these are very inaccurate. A better choice for 
confirming their trajectory character is to verify that they are at a distance 
e of the intersection of the stable and unstable manifolds. The trajectory 
DHTvK remains distinguished from March 5 to May 11, 2003. Outside this 
time interval the path of limit coordinates although still exists for a short 
period, it is not longer a trajectory and eventually the path is lost. This 
kind of behavior coincides with that described in Madrid fc Mancho(2009)] 
for a similar DHT in a highly aperiodic fiow. At the western end there 
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Figure 3: Path of limit coordinates for DHTvk and DHT^; a) longitude vs. 
time; b) latitude vs. time. 

exists a second DHT labelled as DHT,)^, which remains distinguished in an 
approximately complementary time interval, between May 10 and June 1, 
2003. The distinguished property of DHT^ ends similarly to DHTvi/- Panels 
in Fig. |4] show the path of limit coordinates at the East. The trajectory 
DHT^ remains distinguished between March 25 and June 24, 2003. 

The fact that paths of hmit coordinates are close to trajectories, at least 
for a certain interval of time, indicate that the minima of M provide La- 
grangian information. When the motion described in Eqs. (j6])-(j7]) is ex- 
pressed in a new reference frame, trajectories are related to the ones in the 
old frame by the same coordinate transformation. A question then is: do 
trajectories hold in the new frame the distinguished property?. The answer 
is yes. An example is reported for instance in Madrid &: Mancho(2009)] , 



where limit coordinates for the rotating and non rotating Duffing equation 
are related by the appropriate coordinate transformation. 

Once appropriate DHTs have been identified we compute its stable and 
unstable manifolds. A novel technique to compute stable and unstable mani- 
folds of hyperbolic trajectories in aperiodic fiows is developed by Mancho et al.(2003)] . 



In Mancho et al. (2004) , Mancho et al.(2008)| this manifold computation is 



successfully applied to quite realistic fiows. We will now apply these same 
algorithms to compute stable and unstable manifolds of DHT in altimeter 
data sets. In our study we compute the unstable manifold of DHTvf, and 
the stable manifold of DHT^;. The manifold computation requires as an in- 
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Figure 4: Path of limit coordinates for DHT^;; a) longitude vs. time; b) 
latitude vs. time. 

put the position of the DHT at a given time and also the direction of the 
unstable or stable subspace at that time. In this work we obtain this in a 
different way to that proposed in Mancho et al.(2004) . We take advantage 



of the structure of the function M depicted in Fig. |2, which shows strong 
sharp features. In Mendoza &: Mancho(2009)] these lines are reported to be 



an advected structure of Eqs. (l6j)-([7|), thus confirming they are aligned with 
the stable and unstable subspaces. We use these directions as input for the 
unstable or stable subspaces in the manifold algorithm. 

Fig. [5] shows the computed manifolds and DHTs on May 4, 2003. Almost 
every distinguishable line in this figure contains numerous foldings of each 
manifold, thus confirming how intricate they may be. In our computations 
the unstable manifold at a time is made of trajectories that at time to, 
to < were at a small segment aligned with the unstable subspace of the 
DHT. This is a finite time version of the asymptotic condition required for 
unstable manifolds in infinite time systems. Similarly the stable manifold 
at a time t^, is made of trajectories that at time ti, ti > are in a small 
segment aligned with the stable subspace of the DHT. Other approaches such 
as FTLE or FSLE compute manifolds at a given time as ridges of a scalar 
field, providing pieces of curves that are material curves. However links 
between pieces of curves are difficult to establish as they fade away. The 
intricate curves of Fig. |5] are linked curves (due to the asymptotic condition 
imposed in their computation) and this is an advantage over other methods. 
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Figure 5: Stable (blue) manifold of DHT^; at the eastern end of the current 
and unstable (red) manifold of DHT^y at the west. 

Manifolds in Fig. |5] however are too complex to concrete transport problems 
so that in next section we will extract from them pieces with the appropriate 
dynamical information. 

5 The turnstile mechanism 

The turnstile mechanism, extensively used and explained in the literature 



Malhotra and Wiggins(1998) , Rom-Kedar et al.(1990)| , has been found to 



play a role in transport in several oceanographic contexts [Mancho et al.(2008) 



CouUiette and Wiggins(2001)| . Our manifold computations reveal the pres- 



ence of the turnstile mechanism transporting masses of water across the cur- 
rent of Fig. [T]from north to south and viceversa. We study it for a period of 
54 days from April 3 to May 26, 2003. 

First, we determine a time dependent Lagrangian barrier made of pieces 
of manifolds that separates north from south. Manifolds are made of trajec- 
tories, thus as no particle trajectories can cross them, they are barriers to 
transport. Fig. [6] shows barriers at days April 4 and April 17, by depicting a 
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Figure 6: Boundaries at days April 3 2003 and April 17 2003 constructed 
from a (finite length) segment of the unstable manifold of DHT^y and a (finite 
length) segment of the stable manifold of DHT^;. The boundary intersection 
points are denoted respectively by a and b. 

piece of the unstable manifold of DHT^ and a piece of the stable manifold 
of DHT^;. How to choose these portions from intricate curves such as those 
depicted in Fig. [Sp First, we consider that a manifold has two branches 
separated by the DHT, which constitutes a reference point. The selections 
are relatively short length curves either at one or both sides of the DHT. 
For instance for the unstable manifold case, these curves are traced up by 
all the trajectories that expanded from a small segment along the unstable 
subspace of the DHT during a relatively short interval of time. If segments 
along the unstable subspace are expanded for longer time intervals, then 
manifolds become so intricate as those in Fig. |5} In Fig. |6] the boundary 
intersection points are marked with letters a and b. Intersection points sat- 
isfy the property of invariance, which means that if the stable and unstable 
manifolds intersect at a time in a point, then they intersect for all time, and 
the intersection point is then a trajectory. In order to help us to understand 
the time evolution of lobes, we will depict the position of trajectories a and 
b at different times. 

Fig. [7] shows longer pieces of the unstable and stable manifolds at the 
same days selected in Fig. [6] Manifolds intersect forming regions called 
lobes. It is only the fluid that is inside the lobes that can participate in 
turnstile mechanism. Fig. [7] displays two snapshots showing the evolution 
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Figure 7: Turnstile lobes at days April 3, 2003 and April 17, 2003. The 
intersection trajectories a and b are displayed at both days showing their 
evolution from the DHT^y towards the DHT^;. The green area evolves from 
north to south while the magenta area does from the south to the north. 

of lobes from April 2 to April 17. There it is noticed how the lobe which is 
at the north of the barrier at April 2 crosses to the south of the barrier at 
April 17. Similarly the lobe which is at the south in April 2, crosses to the 
north at April 17. Trajectories a and b are depicted showing that they evolve 
approaching to DHTg. The area surrounded by a lobe has been coloured 
with green if near the western DHT, DHT^y, it was at the north side, and 
magenta if it was at the south. The turnstile mechanism transports magenta 
to the north and green to the south. 

Along the time interval under study, several lobes are formed mixing 
waters of north and south. Fig. |8] contains a time sequence showing the evo- 
lution of several lobes created by the intersection of the stable and unstable 
manifolds. A sequence of trajectories a,6, c, ... on the intersection points is 
depicted. These trajectories evolve from west to east and serve as reference 
to track lobe evolution. The lobes between the intersection trajectories a 
and b of [8^) are not depicted in Fig. [8]d) as they are very thin filaments 
stretched along the DHT^;. In the same way the lobes between the intersec- 
tion trajectories c and d of Fig. [8|d) are not depicted in Fig. [8^) as they 
are very thin filaments stretched along the DHTvi/- We notice that in Fig. 
[sji) the western hyperbolic trajectory is labelled as DHT^ since it is a dif- 
ferent distinguished trajectory to DHT^y that stopped as distinguished on 
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Figure 8: Sequence of lobes between April 17, 2003 and May 26, 2003 mixing 
waters from north to south and viceversa at the Kuroshio current. 



May 11. Despite DHTv(/ has lost its distinguished property before we have 
completed our transport description, the manifold computation beyond this 
time still makes sense as it still is a material surface, and is still asymptotic 
to DHTvK- Although the manifold is not asymptotic to DHTj*^, it is correct 
to say that DHT^ marks a Distinguished Trajectory on the manifold within 
certain accuracy e. Any further description on the turnstile mechanism be- 
tween May 11 and 26 needs to be referred to this new DHT. This kind of 
replacement of the DHT of reference along the invariant manifold has been 
discussed in Mancho et al.(2006)] . The piece of stable manifold in [sji) does 
not show more lobes on May 26, and beyond this day the turnstile mechanism 
is disrupted due to topological transitions in the flow structure. 
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The altimeter-derived estimates of the velocity field have expected er- 
rors, so one question is to establish whether the information obtained from 
the analysis of manifolds is robust to changes in the velocity field or not. The 
infiuence of deviations of the velocity fields on the Lagrangian structures has 
been examined in several works |Haller(2002)[[Hernandez-Carrasco et al.(2009) 



Poje et al.(2009)] . The main conclusion of these works is that Lagrangian 



structures are robust to errors, though conditions for the reliability of pre- 
dictions are discussed. In what concerns to lobes and its role on transport 
it is possible to say that these are formed by intersections of stable and 
unstable manifolds and are typically observed in the neighborhood of a Dis- 
tinguished Hyperbolic Trajectory. Despite the lack of rigorous mathematical 
theorems for aperiodic time dependent dynamical systems it is possible to 
foresee a relation between the vicinity of these DHTs and chaotic saddles. 
These latter persist under small perturbations of the vector field, i.e. they 
are structurally stable Wiggins(2003)| . For this reason we also expect that 
the lobe description near the DHTs is robust under small deviations on the 
measured velocity field. 



6 Conclusions 

We have explored the power of dynamical systems ideas for analyzing al- 
timeter data sets. In particular tools for quantifying Lagrangian transport 
have been successfully applied to oceanographic observations. Our study 
analyses cross-frontal transport in altimeter data sets along a fiow located 
in the area of the Kuroshio current. As a novelty with respect to sim- 
ilar studies we have identified relevant DHTs following the methodology 
proposed in a recent paper by [Madrid fc Mancho(2009)|, instead of usmg 



that proposed in |Ide et al.(2002)||Mancho et al.(2004)||Mancho et al.(2006)j 
Mancho et al.(2008)| . We have been able to detect relevant DHTs and we 



confirm that DHTs hold the distinguished property only for finite time in- 
tervals. Identification and computation of DHTs is by itself an important 
subject, since they organize the fiow in the area and, because of this and 
of the sensitivity of the trajectories in their neighborhood, they are candi- 
dates for launch locations in efficiently designed drifter release experiments 
Poje et al (2002)1 [Molcard et al.(2006)] . 



We have computed stable and unstable manifolds of the relevant DHTs. 
A time dependent Lagrangian 'barrier' is defined from pieces of the stable 
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and unstable manifolds and transport across it is described by lobe motion. 
It is found that the turnstile mechanism is at work in this observational flow, 
leading to a filamentous transport near the hyperbolic trajectories. This 
mechanism survives between April 17, 2003 and May 26, 2003. After this 
date transitions in the flow take place, that are not well understood. 

In summary this article explores the use, in a general oceanographic prob- 
lem, of ideas and techniques associated with chaotic transport and dynamical 
systems, trying to shed some light on the mechanisms that lead to messy par- 
ticle trajectories over the ocean. We expect that this tools and information 
will be of interest in the oceanographic community. 
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